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ometry when a non-trivial measure term is included in the gravitational path integral. Within 
our extended region we find a large scale spectral dimension of D s (a — > °°) = 4.04 ± 0.26 and 
a Hausdorff dimension that is consistent with Dh = 4 from finite size scaling. We find that the 
short distance spectral dimension is D s (a — > 0) w 3/2, which may resolve the tension between 
asymptotic safety and holographic entropy scaling. 
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1. Introduction 

The quantization of gravity is one of the great outstanding problems in theoretical physics. A 
straightforward implementation of general relativity as a perturbative quantum field theory is not 
renormalizable by power counting. Although it is still possible to formulate gravity as an effective 
field theory at low energies [[I]], at each order in the perturbative expansion, new divergences appear, 
requiring an infinite number of couplings that encode the physics at high energy scales, leading 
to a loss of predictive power. However, the conjecture that gravity could be nonperturbatively 
renormalizable was made by Weinberg in Ref. [§], where it was termed the asymptotic safety 
scenario. In this scenario, the renormalization group flow of couplings end at a non-trivial fixed 
point in the ultraviolet, with a finite dimensional critical surface. Gravity, therefore, would be safe 
from ultraviolet divergences and describable in terms of a finite number of parameters, making it 
effectively renormalizable when formulated nonperturbatively. 

A lattice formulation allows one to define a gravitational path integral that can be studied 
nonperturbatively. One of the original formulations of lattice gravity is Euclidean dynamical trian- 
gulations (EDT) [||, which defines a spacetime of locally flat n-simplices of fixed edge length, 
where an n-simplex is the n-dimensional analogue of a triangle. An ensemble can be generated 
using the Metropolis algorithm with a set of local update moves. The partition function we use is 



N 2 



e~ s *, (1.1) 



where the product is over all 2-simplices, and 6 (tj) is the order of the 2-simplex7, i.e. the number 
of four-simplices to which the triangle belongs. The term in square brackets is a non-trivial measure 
term. We vary the free parameter /3 as an additional independent coupling constant in the bare 
lattice action [Qj. The Einstein-Regge action is given by, 

S E = -K 1 N 2 + K A N A , (1.2) 

where K2 and K4 are related to the bare Newton's constant and the cosmological constant, respec- 
tively. K4 must be tuned to its critical value such that an infinite volume limit can be taken This 
leaves a 2-dimensional parameter space that can be explored by varying K2 and j3. 

Early studies of EDT typically included only two couplings in the bare lattice action, K2 and 
K4, and they typically did not allow for a non-trivial measure term. The original EDT model ran 
into significant problems. The parameter space of couplings contained just two phases, neither of 
which resembled 4-dimensional semi-classical general relativity, and the two phases were seper- 
ated by a first order critical point, making it unlikely that one could take a continuum limit. In 
response to these problems a causality condition was added by Ambjorn and Loll [Q]. The causal- 
ity condition distinguishes between space-like and time-like links on the lattice so that an explicit 
foliation of the lattice into space-like hypersufaces can be introduced, where all hypersurfaces have 
the same topology. Only geometries that can be foliated in this way are included in the ensemble 
of triangulations that define the path integral measure. A 4-dimensional de Sitter like phase was 
shown to emerge from causal dynamical triangulations (CDT) [0][@]. 

In this work we revisit the original Euclidean theory with a non-trivial measure term for two 
reasons. One is that the CDT restriction to a fixed foliation is potentially at odds with general 
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covariance, though it may amount to a choice of gauge. We would like to see whether the successes 
of the CDT formulation can be reproduced using the explicitly covariant EDT formulation. Second, 
renormalization group studies [KJ suggest that the ultraviolet critical surface is 3 -dimensional, 
providing the original motivation for revisiting EDT with a third parameter in the bare lattice action 



11]. We point out that the CDT action also includes a third parameter, an asymmetry parameter 



between space and time. The results from our EDT approach obtained so far suggest that the 
causality constraint of CDT is not necessary to obtain a good four-dimensional semiclassical limit 
[HI. 



2. Phase Diagram 

The new parameter j8 associated with the non-trivial measure term enlarges the parameter 
space in which to explore the phase diagram. The three different regions of the phase diagram are 
labelled collapsed, branched polymer and extended. 




Figure 1 : A schematic of the phase diagram as a function of K2 and j3 . 



The collapsed region is characterised by a fractal Hausdorff dimension D# — > 00. The branched 
polymer phase has a Hausdorff dimension D# = 2 and a spectral dimension D s = 4/3. These two 
regions are unphysical. However, there exists a region in the parameter space with an extended 4- 
dimensional geometry. The solid line AB dividing the extended and branched polymer phases has 
been extensively studied in the literature at the point /3 = 0. At that point the transition is almost 
certaintly first order [13, 14]. The dashed line CD appears to be much softer than the AB transition 
and appears to be a cross-over seperating the collapsed and extended regions of a single phase. 



3. Numerical Implementation and Code Tests 

A d-dimensional simplicial manifold is constructed by gluing d-simplices together along their 
(d-l)-dimensional faces. To each d-simplex there exists a simplex label and a set of combinato- 
rially unique (d+1) vertex labels. The set of combinatorial triangulations was used in most early 
simulations of EDT. However, the constraint of combinatorial uniqueness can be relaxed to include 
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the larger set of degenerate triangulations in which the neighbours of a given simplex are no longer 
unique [15 J. It has been shown numerically that simulations using degenerate triangulations lead 
to a factor of ~ 10 reduction in finite size effects compared to combinatorial triangulations Jl5|]. 
We have made various checks of our code against the literature using combinatorial triangulations 



[16], as well as for degenerate triangulations Qlip, and good agreement was found. 



4. Exploring the Phase Diagram 

Determining the critical behaviour of the phase diagram is of great importance because the 
existence of a second order critical point allows one to define a continuum limit for the theory. The 
Regge curvature R serves as an order parameter for exploring the phase diagram, with a derivative 
with respect to jc 2 given by the susceptability % R (A^A^). 



N 4 



1, 



(4.1) 



Xr(N 2 ,N 4 ) 




where p 



10arccos(l/4) 
2k 



N 4 , 



(4.2) 



and A 7 ,- is the total number of /-simplices. We have studied the phase transi- 



tion line AB and the transition region CD by calculating the average Regge curvature and curvature 
susceptability for fixed values of j8, whilst varying K%. 




(a) 
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Figure 2: Simulation using 4000 degenerate simplices for j3 = and 4000 and 8000 for J3 = — 1.0, whilst 
varying K%, (a) The curvature susceptability as a function of K-j. (b) The average Regge curvature as a 
function of K%. 



Figure |2| shows a sharp peak in the curvature susceptability for j8 = 0, and a less sharp peak for 
j8 = — 1 .0. The behaviour at j3 = — 1 .0 is consistent with a cross-over as the peak in the curvature 
susceptability does not appear to change with volume. 1 

' The transition could be third or higher order, since it is very difficult to distinguish an analytic cross-over from a higher order 
transition using numerical methods. 
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It is well established in the literature [13] that, at the point j3 = 0, along the critical line 



AB, there is a first order phase transition for combinatorial triangulations, and Ref. Q17| ] provides 
evidence that this is also true for degenerate triangulations. We show the Monte Carlo time history 
of No and the corresponding histogram for the ensemble at the peak in %R at K"2 = 1-672 and j3 = 0. 
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Figure 3: Simulation using 4000 degenerate simplices at K2 = 1 .672 and j3 = 0. (a) The number of vertices, 
A^o, as a function of Monte Carlo time, (b) Histogram of A^o values. 

Figure [3a| shows evidence for discontinuous fluctuations between the two metastable values of No 
that are characteristic of a discontinuous phase transition, and Figure [3b| shows the histogram of 
No- A first order phase transition has a signature double Gaussian peak in the histogram of No, as 
seen in Figure RR 



5. The Extended Region 

Simulations investigating finite-size scaling support the fact that there exists a region in pa- 
rameter space where our EDT model has a 4-dimensional extended geometry. We study the finite 



volume scaling behaviour of the three-volume correlator similar to the one introduced in Ref. [ |18| ] 
to study the scaling of CDT, 



Cn 4 (8) = £ 



i, (Nf ce {z)Nf cs {T + 8)) 



(5.1) 



A/4 (t) is the total number of 4-simplices in a spherical shell a geodesic distance % from a ran- 
domly chosen simplex. A/4 is the total number of 4-simplices and the normalization of the cor- 



relator is chosen such that T, s=q Cn 4 (8) = 1. If we rescale d and (8), defining x = 8/N A 
then the universal distribution cn 4 (x) should be independent of the lattice volume, where cw 4 (x) 



N 



Cn 4 S/N t 



One can determine the fractal Hausdorff dimension, D#, as the value that 



leaves ca? 4 (x) invariant under a change in four-volume A/4. Good agreement between c^ 4 (x) at dif- 
ferent volumes A/4 occurs when D# = 4; thus the geometries have Hausdorff dimension close to 4. 
This is illustrated in Figure ^3. 

The spectral dimension, D s , is related to the probability of return, P r (o), for a random walk 
over the ensemble of triangulations after a diffusion steps. It is defined by, 
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Figure 4: (a) Three-volume correlator using 8000, 16000 and 32000 four simplices. (b) The running spectral 
dimension as a function of diffusion time, including a fit function of the form D s — a 



/logfi-(q) 

dlog a 



(5.2) 



Figure [4b] shows a plot of D s as a function of a, including a fit to the function D s = a — sug- 
gested in Ref. Jl9|]. Our preferred fit to the spectral dimension gives D s (a — > °°) = 4.04 ± 0.26 
and D s (a — >0) = 1.457 ± 0.064, where the error includes the statistical error and a systematic 
error associated with varying the fitting range and the fit function added to the statistical error in 
quadrature. Variations in the fit function assume that the spectral dimension approaches a constant 
in the a — > °° limit, and that D s (a) is monotonic. Although this result for the spectral dimension 
contains an estimate of the systematic error, to bring this error fully under control we will need 
to explore finite size effects for D s (a — > °°) and discretization effects for D s (c — > 0). Work in 
progress using ensembles with larger volumes and longer diffusion times still give D s (a — > °°) 
consistent with four. It is especially important to quantify a systematic error associated with dis- 
cretization effects, because our result for D s (a — > 0) is inconsistent with the renormalization group 
result D s (a — > 0) = 2 exactly p0|]. 2 So far, calculations at finer lattice spacings extrapolate to val- 
ues close to D s (a — > 0) = 3/2 and show no sign of increasing or decreasing. 3 



6. Conclusions 



We have presented evidence that EDT with a non-trivial measure term has a region with 4- 
dimensional geometry. This geometry has a large scale spectral dimension of approximately four, 

2 Ref. [|]] argues using renormalization group methods that there will be a long plateau at ~ 4/3 in D s (a) for small values of a 
before D s increases to 2 as a — > 0. Thus, they propose that our value D s (a — > 0) ~ 3/2 may be the result of our using insufficiently 
fine lattice spacings and a non-monotonic behaviour of D s (<j) for small o". 

The relative lattice spacing is determined from a comparison of the running spectral dimension at different values of the bare 
parameters. 
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and a Hausdorff dimension that is consistent with Dh = 4, as determined from finite size scaling. 
These results suggest that the explicitly covariant EDT approach to quantum gravity may be in the 
same universality class as CDT. The small scale spectral dimension of D s (cr — >■ 0) = 1.457 ±0.064 
is consistent with 3/2, and the CDT result D s (a — > 0) = 1.80 ±0.25 [ |l8| ] is also consistent with 



3/2, given the current error estimate. As discussed in our earlier paper Q12|], this result may resolve 
the tension between asymptotic safety and holographic entropy scaling. 

This work was funded by STFC and the Scottish Universities Physics Alliance. Computing 
was done on Scotgrid and the DiRAC facility jointly funded by STFC and BIS. 
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